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ABSTRACT 


The evolution of forces acting on horizontal cylinders Subjected To 
impact by an harmonically oscillating free surface has been investigated 
both theoretically and experimentally. 

| The experiments were carried out in a large U-shaped water tunnel, 
with cylinder diameters ranging from 3.0 Fe 6.5 inches. The results 
have been expressed in terms of two force coefficients. The first is 
the slamming coefficient which expresses the normalized force acting on 
the cylinder at the time of impact; and the second is the maximum drag 
coefficient which occurs when the cylinder is immersed approximately 
1.7/5 diameters in water. The slamming coefficient was found to be equal 
to 7, as predicted theoretically. It was also found that the slamming 
coefficient may be amplified to a value as high as 2m through the 
semirce response of the cylinder. The maximum of the normalized force 
In the drag dominated region was found to be approximately equal to 2.0. 
it decreased with increasing Froude numbers to a value of about unity. 

It 1S recommended that further experiments be conducted with sand- 
roughened cylinders and the dependence of the rate of rise of the impact 


force on the flow parameters be investigated. 





lel. 


EY. 


Nes 


TASe CE COMMENTS 


Wiicememiepeeeen= = - = =-- =----2.+22------- 
eee MT nySiS = = —- —-— — — - - — - Se ge 
EXPERIMENTAL EQUIPMENT AND PROCEDURES - --------- 
Ay RAINE © ce he 
PREIETRBERS = = - - -=---2----------.-.-- 
C. FORCE MEASUREMENTS -- ---------------- 
RCUESMENMOr FLUID VELOCITY - ---=-=-----+.-- 
ERMEFPESSSNING EQUIPMENT - - --=-----+-+------- 
Rea ENMeROCEDURE = - - - - - - =e ---..-- 
PRPS Os DATA - = --------+-+-------+ 
P@ISseIMOr RESULTS = -=2 = --+-=+-----+-+---- 


CONCLUSIONS - -------+-+---+-+---+--+----- 


LIST OF REFERENCES - - - ----+--+--+-+---+-4-4--4-424-- 


INITIAL DISTRIBUTION LIST - ---+-+--+-+-+-------4--°- 


23 
ZS 
eS) 
Se 
Sz 
S72 
33 


34 


US 
“6 


U7 





Figure Page 
|. Geometry for theoretical analysis of wave impact - - ---- - 14 
2. om/dz as a function of 2/D -------------+------ LG 
3. Theoretical values of C_ as a function of 2/D 

for several values of no/A -------+----------- 18 
4. C(t) as a function of T------------ - ee ee ee 20 
5. k Gree as a function of + for t_ = 0.0001 second - - - -~ - 93 
oe «Kk nen e. sear cuncrtion: ot tf fOr We = 0.010 second - - --- - gk 
7. k x(t)/C> as a function of t+ for t. = 0.0195 second - - - - - 25 
en = K x(t) /Ce aoe Tunction of fF for i = 0.023 second - ---- - 26 
oe CK x(4)/C2 as a function of ft. --------------- 98 

lO. Schematic diagram of the U tunnel - - - ----------- 30 

l!. Data trace showing filtered and unfiltered traces - ----- 35 

12, Data trace for a run with double first peak - -------- 36 

mm. Odata trace for a run with decreased first peak - - - ----- ey 

14. Ce (Theoretical), C. (experimental), and buoyancy 

aS a function of 2/D--------+-+--7 7-7-7 - - - ee ee 39 

. i ve Samar runechion Of g0/u, = |/Fr------+-+-+-+-+-+-- 4 

l6. C, as a function of Fr-------------+--+--+-+--- men 


Pio Or a CURES 


h 





NOMENCLATURE 


Amplitude of the wave motion 
Immersed cylinder area 

Normalized immersed cylinder area 
Buoyancy corrected drag coefficient 
Slamming coefficient 

Initial value of the slamming coefficient 
Cylinder diameter 

Instantaneous total force 

Froude number 

Natural frequency (Hz), Qa TA - 
Gravitational acceleration 
Keulegan-Carpenter number 

Length of cylinder 

Mass 

Added mass 

Normalized added mass 

Pressure 

Reynolds number 

Radius of cylinder 

Distance between pressure taps 
Wave period 

Time 


Rise time 





Maximum free surface velocity 
Instantaneous velocity 

Initial velocity 

Depth of cylinder immersion 

Normalized depth of cylinder immersion 
Damping coefficient 


Instantaneous height of wave surface above 
mean water level 


Height of cylinder with respect to mean water level 
Fluid density 

Velocity potential 

Natural frequency (rad/sec) 


Viscosity 
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» INTRODUCTION 


Information concerning the forces acting on bluff bodies subjected 
to wave slamming is of significant importance in Mechanical Engineering, 
Ocean Engineering, and Naval Architecture. The design of structures 
which must survive In a wave environment is dependent on a knowledge of 
the forces which occur at impact, as wel! as on the dynamic response of 
The system. Two typical examples include the structural! members of off- 
Shore drilling platforms at the splash zone and the often-encountered 
Slamming motion of ships. 

The general problem of hydrodynamic impact has been studied exten- 
sively L!] motivated in part, no doubt, by its importance in ordnance 
and missile technology. Extensive mathematica! models have been developed 
for cases of simple geometry such as spheres and wedges,and these models 
have been well supported by experiment. Unfortunately, the special case 
of wave impact has not been studied as extensively. Kaplan and Silbert 
L2] developed a solution for the forces acting on a cylinder from the in- 
stant of impact to full immersion. Dalton and Nash [3] conducted slam- 
ming experiments with a 0.5 inch diameter cylinder with smal! amp!itude 
waves generated in a laboratory tank. Their data exhibited large scatter 
and showed no particular correlation with either the predictions of the 
hydrodynamic theory or identifiable wave parameters. |t appears that the 
evaluation of the slamming effects with wavy flows is extremely difficult 
partly because of the !imited range of wave amplitudes that can be achieved 
and partly because of the difficulty of measuring the fluid velocities 


Seine instant of impact. 
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In essence, the present research was undertaken to extend the work of 
Dalton and Nash and to improve the experimental method by using a large 
U-shaped water tunnel. In addition, it was intended to demonstrate a cor- 
relation between theory and experiment by taking into account the dynamic 
response of the system. Specifically, the goals of the project were as 
fol lows: 

A. to examine the existing theoretical models in determining wave 

Slam forces on circular cylinders; 

B. to furnish data, obtained under controlled laboratory conditions, 
about the forces acting on smooth circular cylinders subjected to 
impact with an harmonically oscillating water surface; 

C. to determine the relative importance of the inertia dominated and 
drag dominated forces during wave impact; and 

D. to correlate such data with respect to identifiable wave parameters 
such as the Froude number, Fr; the Keulegan-Carpenter number, K; 


and the Reynolds number, Re. 
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eee eeOneinicat ANNE YSIS 


The traditional approach used in the design of off-shore structures 
involves the classical Morison equation to determine the forces due to 
wave motion. However, wave impact is generally defined as the early 
stages of the penetration of a solid body into a wave surface. At the 
instant of contact the fluid in the vicinity of the body undergoes large 
accelerations which give rise to large forces. As the body becomes more 
fully immersed, forces due to viscous drag and separation effects become 
predominant. Thus the inertia and drag coefficients used in the Morison 
equation are not constant,and the problem becomes very difficult, even 
for simple geometries. 

The general case of hydrodynamic impact is usually described by using 
incompressible potential flow theory. For the case of a moving body with 
mass,M,and velocity,v., impacting a quiescent surface, the system momentum 
is ey Neglecting non-conservative forces, the momentum of the system 
iS unchanged during penetration. However, the mass of the system 
Beemeases due To the fluid which is set in motion in The vicinity of the 
body. Also known as "added mass", m, it results in reduction of the 
velocity. Thus the system momentum after penetration is (M+m) v = Mv oe 
The impact force at any instant is a function of m and dm/dt. The 
solution, therefore, requires knowledge of the added mass and its time 
derivative. 

it should be noted that viscous effects, being non-conservative, wil 


alter this somewhat. Additionally, high speed entry (as witn projectiles) 
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may require that compressibility effects of both the fluid and the air 
above it be taken into account. However, for low speed entry, these 
effects are generally neglected,and the fluid motion is described by a 
velocity potential, , which satisfies y od = 0. The added mass term 
can then be calculated. The problem remains extremely difficult, though, 
because the free surface tends to "pile up™ around the body. This 
results in an unsteady flow problem with a time dependent boundary condi- 
tion (the free surface). Approximate solutions have been obtained [4] 
for simple shapes (spheres, cones, cylinders, wedges, flat plates, etc.) 
by assuming a constant plane free surface and then applying a correction 
mor the distorrion. 

Kaplan and Silbert [2] developed the following relation for the force 
acting on a horizontal circular cylinder by a wave system which propa- 


Manes normal to it, (See Fig. |.): 


F -- om eZ 
i= Og A, sae Got ae aia @ Aj) n ae n ole. 


in which F represents the force acting on the cylinder; L, the length 

of the cylinder; 9, the fluid density; g, the gravitational acceleration; 
A. the immersed area; m, the added mass per unit cylinder length; n, The 
instantaneous height of the wave surface above the mean water level; and 
Z, the instantaneous depth of cylinder immerston. The first and second 
derivatives of n with respect to the time are denoted by fi and nj. The 


added mass is given by Taylor [5] as: 


5 
7 er is (|=cos G oe 


aC Tae a z(1-cos See Sin Go =o) | C25 


in which r represents the radius of the cylinder, and 8 is defined as 


ewewn in Fig. |. 


Ais) 





| Aj 
7 
M.W.L. 


Fig. |. Geometry for theoretical analysis Of wave impact. 
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The motion of the free surface, n, is related to the maximum ampli- 


tude by: 


a 2 A sie eee (3) 


where A and T represent the amplitude and period of the free surface 


oscillations. Equation (|) can also be written in the form of a slamming 


Boetticient, Ce =n Oe ue as: 


Pe ee gy Dig So ee 2M cost Sh + (4) 
S | 2 i A sl = ae 
U Oz 
m 
where 
A. = A./r? , m= m/pr 
yaw, = 2 A/T 


With z = r(l=cos $), Equation (2) can be used to relate the added 
mass to the depth of immersion. Figure 2 is a plot of dm/dz as a func- 
tion of z/D. Cleariy, 3m/az begins with an initial value of 7m and drops 
mapidly. the quantity gee is related to the Froude number, Fr, by 


= U */2gr. ioieee Cn tGty 5 iy 


A 
dZ 
The rate of change of the normalized added mass with z depends on 6 


and hence on the time measured from the instant of impact (See Fig. 2). 


For example, for the case where ae = Orang t= 0, om/oz = 1 and Ce 


T 
=3 


Consequently, the for the particular case under consideration, os at the 
instant of impact does not depend on either the size of the cylinder or 


the flow parameters. 


ils 
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The fact that the cylindrical members of a structure in the splash 
zone are not necessarily located just at the mean water level requires 
the determination of the particular value of Ms for which the slamming 
force at the time of impact is a maximum. I+ can be demonstrated through 
the use of Equation (4) that the maximum impact force occurs for the 
case of i O. For this purpose Equation (4) was evaluated with the aid 
of a computer program which allowed variations in r/A and no/A- C. was 
then plotted as a function of z/D from zero to unity. Figure 3 is an 
example of one of these plots and is typical for most cases. The variables 


neTe; 


mA = 0.125 


n/A Cee. 0245 and 0.3 
As can be seen readily from Fig. 3, C. is largest at z/D = 0 for n/A =0.0 
and starts at a value of m7 and drops to a minimum at z/D of approximately 
OE). it rises again as z/D approaches |.0. Thus Equation (4) indi- 
cates that Co, and consequently the impact force, is of an impulsive 
nature beginning with a finite value (not zero) at the instant of impact. 
Since viscous effects are neglected, one would expect the solution to 
deviate from the actual situation as the cylinder becomes more fully 
immersed. Where this becomes the case can only be determined by experiment. 
However, one can assume that the solution would be valid for smal | 
depths of immersion, which is the case for the instant of impact up to 
meat 0.02 or so. With this restriction, A. and m are small, as is 
Sin 2mt/T. Additionally, wage 2mt/T is very nearly equal to 1.0. There- 
fore C. reduces to: 


om -_ —_= 
as Cam/9Z)4_4 


iy; 





Theoretical values of C. a 
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It is not realistic to assume that the impact force rises from zero 
to wm instantaneously. Several factors, specifically the compressibility 
of the air between the cylinder and water surface, entrapped gases in 
the water, and surface irregularities would account for some finite rise 
time. Nonetheless, the rise time can be expected to be short, i.e. in 
the order of milliseconds. The exact nature of the rise is an interest- 
ing question for further study. However, in this analysis Ce IS assumed 
To vary linearly during a rise period, Tee POS 1a es Ce drops from 7 
as Om/dz. Figure 4 is a representation of this assumption. Exactly how 
long of a time interval e is will be discussed later. 

The realization that the impact force is of an impulsive nature requires 
consideration of the fact that this force does not act on a perfectly rigid 
body, but rather on a cylinder which is supported elastically. The re- 
sponse of such a system approaches that of a rigid body only if its 
natural frequency approaches infinity. Additionally, the response of The 
system to an impulsive force is heavily dependent on the exact nature of 
The force itself as well as on the system natural frequency. 

In general, the instantaneous displacement, x(t),of a system of mass, 
M,with a spring constant,k,is given by (See e.g. L6]): 

a 
x(t) -f mie @ Chae )de (6) 


O 
in which € represents a dummy variable; F(t), the driving force; and 

g(t), the response to a unit step excitation. Equation (6) is readily 
recognized as the DuHame!l superposition integral, which can be expanded 
as follows for g(t) = 


| e J 
——- S/N W T: 
Ma n 
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os | , 
x(7) = Fas °° ot f Seem siig w tat 
O 
a 
| 
+ 
My $1" ut f Fon cos Ww Tdt FO) 
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in which W is the natural frequency of the cylinder and supports. 
Letting F(t) = Flos and changing the variable of integration to 


9 = wit, one obtains 


9/ 8/u) 


KAY = = cos a ft sin 6d8 + sin 6 f te 00s ade (8) 


‘e) 
oO 


. ; 2 ; 
In which k = W Pane Spring consnani. Wme férm on the left-hand side 


can be interpreted as the ratio of the instantaneous force sensed by the 


Supports of The cylinder to the actual mean force acting on the cylinder. 


In other words, the response of the system at Time + is k x(t) which may 


be different from Foe For the simple case in which F(T) is a step 


function as shown below, 


equation (8) reduces to: 


ax (1 = cos wt) (9) 


oO 





It is apparent from Equation (9) that the instantaneous force sensed by 
the system can be anywhere from zero to two times the actual mean force. 

lf one assumes, as before, that the impact force is as described by 
Fig. 4, then Equation (8) must be evaluated using C.(1) Sai Gd) acid = 


tionally, damping can be taken into account by writing g(t) as [6]: 
3h) = e S45" sin wt (10) 


where [ is the damping coefficient. Thus Equation (8) can be rewritten 


by replacing the forcing function, F(t), by C (8) as: 


8/u 8/0) 
kx) 2 cose f (a)e —"singde + cose f C (a)e o%cos6dé Gia) 
oo : : 
S O O 


Equation (II) was solved with the aid of a simple computer program which 


used a trapezoidal integration of variable step. Values of w were 


taken as 358 sec. and 628 sec! which corresponded to the measured values 


a7 A 


of 0 for a 6.0 inch and a 3.0 inch diameter cylinder respectively, as 
discussed later. The damping coefficient was 0.014 in both cases, also 
corresponding to measured values. The rise time was varied from zero To 


approximately 0.025 seconds. Figures 5 through 8 are representative plots 
f(T) 


Cue 
Figure > represents Equation (||) plotted for a rise time of Q.000I 


of LG ticmomlmen cylinder, tom various values of Th. 


maeema. Fig. 6, for a rise Time of 0.0100 second; Fig. /, for a rise 


Time of 0.0195 second; and Fig. 8, for a rise time of 0.0230 second. 
K CP) 
O 


& 
of approximately |.7 at t = 0.001 seconds and then drops off. As the rise 


Note that for a very short rise time the value of reaches a value 


Time increases, the value of the first peak drops off, with a "double 


peak" appearing at a rise time of 0.0195 seconds. 
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t. = 0.010 second 


dal ecaye (3) kx(t)/Co aoeaehunectwam of t 


for ea —'O20 POs Second.. 


24 





1S ce = 0.0195 second 





0 On02 0.04 0.06 


|S teas al kx(t)/C° as a function of t 


for cae 0.0195 second. 
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i ees oe = 0.023 second 





19 ot fees kx (tye? as a function of t 


for Gs = 0.023 second. 
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The interpretation of the above results indicates that, depending on 
the rise time of C.(T), values of the apparent slamming coefficient, om, 
may !ie between roughly 0.5 and 1.5 of the theoretical instantaneous 
value of 7. Again this applies only to a smal! depth of immersion, cor- 
responding to the initial moments of impact. The significance of this 
pae7 is that values of Co measured in the laboratory may show wide scat- 
ter depending on the rise time for any given experiment. Additionally, 
if the surface is not perfectly plane, rise times may vary from experiment 
to experiment resulting in an apparent non-repeatability. 

mireure’ 9 is a plot of k Gey ee as a function of a Te for both the 
first and the second peak, using w Vendes: TO bOr a 6.0 iimch and a 5,0 
inch cylinder. From this figure, the drop in k Gare. for the first 
peak with increasing rise time is apparent. It is also evident that the 
effect of We is confined to a narrow range of ie ues values smaller than 
about 0.5. The force acting on the cylinder is amplified by the dynamic 
response of the system for a value of i ih smaller than 0.9 and 
attenuated for i ie (Qe oelceagceatiaieOes..) Ine SeCOnd peak Occurs 
emily tor a ie values larger than !.09. The value of k SCPC. increases 


from 0.84 to about unity as se ie increases from 1.09 to |.5. 
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lll. EXPERIMENTAL EQUIPMENT AND PROCEDURES 


The equipment used to generate the harmonically oscillating flow has 
been used extensively at this facility over the past two years. The 
apparatus is described in reference [7]. The salient features, 
as well as the adaptation for this work, are briefly described in the 


following. 


fee TUNNEL 

Figure 10 is a schematic diagram of the U tunnel. A butterfly valve 
arrangement at the top of the left leg allows that side to be completely 
sealed off. Compressed air can then be admitted until the water level 
is driven to the desired height in The opposite leg. Rapid opening of 
the butterfly valve allows the water to oscillate at a natural frequency 
determined by the geometry of the tunnel. The tunnel contains approxi- 
mately 5000 gallons of water and the period, T, iS 5.5 seconds. Ampli- 
tudes up to 4.0 ft can be generated. 

Previous work with this Tunnel has demonstrated its capability in 
Studies involving oscillating flow. Measurements have confirmed that 
Separation does not occur at the corners and there are no discontinuities 
where the individual sections are joined, which might otherwise lead To 
disturbances in the flow. These factors were important to this project 


Since uniformity of the free surface was a major concern. 


Eee ves) CYLINDERS 
Three, five, six, and six and one half inch diameter aluminum 


cylinders were used in measuring impact forces. Measuring 2.98 ft in 
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length, they were constructed of aluminum pipe, turned to the final 
diameter, and polished to assure a hydrodynamically smooth surface. 
Caps fitted to each end prevented water leakage into the cylinder. 
Each cap contained a double precision ball-bearing mounted flush with 
the face. As will be noted later, the force transducers were attached 
to the cylinder via these bearings,and the purpose was to allow free 


rotation of the cylinder. 


SeeeerOrceE MEASUREMENTS 

Two cantilever-beam force transducers were used to measure the in- 
stantaneous in-line (drag) and transverse (lift) forces. The output 
consisted of an electrical signal from strain gages attached to the 
beam. The rounded end of the transducer assembly was designed to fit 
snugly in the bearings in the ends of the cylinders. Neoprene "0" 
rings were added to the end of the transducers to provide axial position- 
ing of the cylinder within the tunnel. 

Calibration of the transducers was accomplished by hanging a load 
from the center of the cylinder. This not only established the level 
of the electrical output for a known load but also reaffirmed that the 
Transducers were correctly oriented. 

As with other elements of the experimental apparatus, these tTrans- 
ducers had been used frequently at this facility for over two years and 
the description of their construction may be found in reference [7]. 
No changes have been noted in the calibration of these transducers since 
Their installation, a fact which establishes some degree of confidence 
in their use. 

Initial experimental efforts established that the transverse forces 
were very small compared to the in-line forces,and consequently measure- 


ment of these forces was discontinued early in the study. 
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D. MEASUREMENT OF FLUID VELOCITY 

Two pressure taps, one in each leg of the tunnel, were connected by 
tubing to a differential pressure transducer located midway between the 
taps. The acceleration of the fluid can be determined from Ap = ps a 
in which Ap represents the pressure differential; po, the fluid density; 
Ss, the distance between the Taps; and — The instantaneous acceleration 


of the fluid. The maximum velocity and amplitude are related to the 


acceleration by: 


(=) =—U = (=) A (12) 


E. RECORDING EQUIPMENT 

A three channel strip chart recorder was used to record simu|ltaneous} y 
the output of the in-line force transducer and the differential pressure 
transducer. The third channel of the recorder was also used to record 
the force measurements, but the signal was first passed through a low- 
pass filter to remove the oscillation of the cylinder at its natural fre- 
Guency, thus providing an "average" force record. 

Inasmuch as the impact force to be measured consisted of a rather 
Short rise time, the bandwidth of the recorder was of some importance. 
The frequency response of this recorder was flat to 100 Hz. In order to 
verify the response under the test conditions, several experimental runs 
were made in which the force was both recorded using the strip chart re- 
corder and a tape recorder simultaneously. The taped signal was then 
fed back to the strip chart recorder at |/4 speed. This method assured 


that the frequency response of the strip chart recorder was adequate. 
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F. EXPERIMENTAL PROCEDURE 

The success or failure of any experimental work depends on careful 
attention to those factors which may have a significant effect on the 
measurements being taken. For this reason,extreme care was taken with 
all flow parameters. 

Prior TO filling the Tunnel, the cylinder was struck lightly,and the 
natural frequency and damping coefficient were determined. The Tunnel 
was then filled, immersing the cylinder several diameters, and again the 
natural frequency was recorded. Thus the value of w was obtained both 
in air and in water. In the calculations the value of Oo, in air was 
used since at impact and shortly thereafter, the cylinder behaves as 
though it were not immersed. 

Since the height of the mean water level with respect to the bottom 
of the cylinder was shown by the analysis to be important, the water 
level in the tunnel was adjusted carefully by slowly filling or emptying 
The tunnel until a slight ripple occurred at the water surface due to 
contact withsthe cylinder. This method assured that the mean water level 
coincided with the cylinder bottom surface. 

Of paramount importance was the condition of the free surface at 
The instant of impact with the cylinder. Ideally the surface should be 
perfectly plane with no surface disturbances. It is this requirement 
That lends a certain degree of difficulty in the experimental procedure. 
Until the fluid in the tunnel has completed one cycle or so, the period 
of oscillation is not established. This requires that the cylinder be 
immersed at least once before the period of oscillation has settled down, 
resulting in some surface disturbance. The observations of the wave 


impact indicated that one slam generated minimal surface disturbances, 
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with the disturbances increasing with subsequent impacts. For this 
reason, all data were recorded for the first impact only,and successive 
impacts were discounted. Approximately ten minutes were allowed between 
each experimental run in order to allow the water surface to reach a 
quiescent state. 

The stripchart recorder was run at a speed of 200 mm per second, each 
division representing 0.005 second. Figure ||] is a typical example of 
The data obtained. The top trace is the filtered force trace,and the 
bottom one is the unfiltered force trace. Figures 12 and |I3 show addi- 


mien aieeuatt|lierea force fraces for various initial flow conditions. 


eee eOUCTION OF DATA 


Two force coefficients were defined in this tnvestigation. They are: 


Ze 


Ce = ee Cl? 
pDLU, 
and 
C= A - EE (14) 
oDLU 2U 
m m 


The first of these refers to the slamming coefficient as defined by 
Equation (4). The second coefficient represents the second maximum of 
The normalized force with the buoyancy subtracted. This maximum occurs 
after the cylinder is fully immersed. 

As mentioned previously, the stripchart recorder was run at 200 mm 
per second, thus each mm corresponded to 0.005 second. The measured 
force was then read each 0.005 second up to 0.5 seconds or so, depending 


On the diameter of the cylinder. These points were then evaluated with 
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Data trace for a run with decreased first peak. 


Pages: 





the aid of a computer program which solved Equation (4) and plotted 

oS as a function of z/D. Figure 14 is typical of such plots, which shows 
the measured value of Co, The theoretical value of Ce from Equation (4), 
and the normalized buoyant force. 

The filtered force trace obscured the magnitude of the initial impact. 
However, after the initial peak, the filtered trace was found to corre- 
spond to the average value of the unfiltered trace. Therefore, the data 
points during the early stage of impact were read from the unfiltered 
trace,and the remainder were read from the filtered trace. The main rea- 
son for this is the obvious difficulty in interpreting the average value 


which iS obscured by the oscillation of the cylinder. 
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vem ON SGUSo ION OF RESULTS 


The initial value of the slamming coefficient, ost was eas as a 
U 


function of various parameters such as the Froude number, Fr = il ; The 
Reynolds number, Re = U_D/v ; and the Keulegan-Carpenter number, K = 
U_T/D. These plots have revealed that there is no identifiable corre- 
lation between C. and the said parameters. In fact, it is the realiza- 
tion of this lack of correlation that led to the consideration of the rise 
time, le Subsequently, the amplification or attenuation of C. was calcu- 
lated by dividing ca by its theoretical value of T Then the correspond- 
ing i De values were determined from Fig. 9. Figure 15 is a plot of 

ia ge as a function of the inverse Froude number g0/u_*. Even though 
there is considerable scatter, this figure shows that for small values 

of Uy or large vaiues of D, the rise time is larger and thus yields 
smaller amplification (See Fig. 9,)., For larger impact velocities or 

for smaller cylinders, the rise time is considerably smaller, and +he 
amplification factor approaches |./7. |+ should be noted 

once more That surface irregularities, dissolved air, and other experi- 
mental uncertainties do not permit a more accurate determination of the 
rise Time. One must add that the rise time may vary also with the 
physical characteristics of the fluid above the free surface, the angle 
of inclination of the cylinder, the length to diameter ratio, etc. Be 
That as it may, the Froude number appears to be the dominant parameter 

as in all free surface flow problems. 


The physical reasoning behind the correlation of the normalized rise 


time, f. t_»with gb/u_ is as follows. The disturbances in the free 
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surface during the first passage of the flow increase with increasing 
diameter. On the other hand, the larger the amplitude, the longer the 
time for the disturbances to attenuate. Thus, the rise time decreases 
with increasing amplitude. Finally, the larger the acceleration of the 
mow i.e. A(2n/T)* relative to the gravitational acceleration, the 
faster is the rate of attenuation of the disturbances [8,9]. A simple 
dimensional analysis shows, therefore, that Lie Tis. increases with 
gD/A(T“/A). Since U_ is proportional to A/T, one finds that f ihe 

a function of g0/u*. 

A careful comparison of the measured and calculated force traces 
during the early stages of impact lends ample credence to the idea of 
rise time. Figures 5 through 8 and I! through |I3 show that for very 
Short rise times there is a single sharp peak In The force trace. For 
intermediate rise times (See Figures 7 and 12.) a double peak occurs. 
Finally, for larger rise times (See Figures 8 and 13.) The first peak 
practically disappears, and the response to the slamming force is 
attenuated by the system. 

It ts event from the foregoing that the determination of The initial 
Slamming force is extremely difficult because of its magnification or 
attenuation by the dynamic response of the system. The evidence presented 
above shows That the normalized slamming force should be assumed equal To 
Tm and that the response of the system be determined through a straight- 
forward vibration analysis. Such an analysis will yield the largest 
force amplification for a normalized rise time of i a Ce ae ae) 6 
design purposes, the maximum slamming coefficient may be as large as 1|./T. 

Following the initial impact, the net force acting on the cylinder 


begins to decrease since the dm/3z contribution decreases. During 
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this stage of the flow the cylinder undergoes damped oscillations at its 
natural frequency. As the free surface rises further, the buoyant force 
increases,and the separation effects begin to give rise to larger drag 
forces. The buoyancy subtracted fluid force reaches its maximum at z/D 
Values from about |.5 to 2.0. Even though the present flow situation, 
in which there is a free surface,cannot be directly compared with an 
impulsive flow about a circular cylinder, the rise of the drag force to 
@ maximum at z/D — 1.75 is very much like the rise of the drag coeffi- 
cient to a maximum at a relative fluid displacement of 2.0 in impulsive 
flow [|]O]. The said rise in the drag coefficient is because of the for- 
mation of two symmetrical vortices behind the cylinder. As the motion 
Gondinues,the vortices become asymmetrical and shed alternatingly. 

The maximum force corresponding to The second peak has been 
corrected for buoyancy, normalized by = pU_,DL,and plotted as a 
function of the Froude number in Figure 16. Two facts are apparent: 
First, the scatter in the data is considerably less than that correspond- 
ing to the initial impact; second, C, [See Equation (14).] decreases 
with increasing Froude number to a value of about unity and remains 


Geaiewy COnNStant for a!l Froude numbers from 0.5 to 2.0. 
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Vo CONCEUS IONS 


The theoretical and experimental investigation of flow impact on 


circular cylinders warranted the following conclusions; 


The dynamic response of the system is as important as the impact 
force and one cannot be determined without taking the other into 
consideration. 

The initial value of the slamming coefficient is essentially equal 
to its Theoretical value of 7, and the system response may be 
amplified or attenuated depending upon its dynamic characteristic. 
The normalized rise Time appears to be a function of the Froude 
number. 

Following the impact, the cylinder undergoes damped oscillations 
at its natural frequency, 

The buoyancy corrected normalized force in the drag dominated 
region reaches a maximum at a relative fluid displacement of 
moun 2/D — e798 subsequently, the shedding of vortices as 

well as the deceleration of flow reduces the maximum drag coefTfi- 
cient to below unity. 

Since it does not appear that the rise time can be determined 
Theoretically, the dynamic response of the system should be 
analyzed using the theoretical value of the slamming coefficient 


and the experimentally determined rise times. 
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